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Completeness is Unnecessary for Fast Nonlinear Quantum Search 
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Although strongly regular graphs and the hypercube are not complete, they are “sufficiently 
complete” such that a randomly walking quantum particle asymptotically searches on them in the 
same Q(\fN) time as on the complete graph, the latter of which is precisely Grover’s algorithm. We 
show that physically realistic nonlinearities of the form /(|V’| 2 )V' can speed up search on sufficiently 
complete graphs, depending on the nonlinearity and graph. Thus nonlinear (quantum) computation 
can retain its power even when a degree of noncompleteness is introduced. 
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I. INTRODUCTION 

Although Grover’s quantum search algorithm jlj was 
originally proposed as a digital algorithm, where the state 
of the system evolves in discrete-time when acted upon 
by quantum gates, it can also be formulated as an ana¬ 
log algorithm, where the system evolves in continuous¬ 
time by Schrodinger’s equation with some Hamiltonian. 
This continuous-time analogue of Grover’s algorithm was 
first proposed by Farhi and Gutmann [2] , and Childs and 
Goldstone [3] later provided a physical, intuitive way to 
interpret it as a quantum particle randomly walking on 
a complete graph with N vertices. Labeling the vertices 
{|0>, 11),..., \N —1)}, the system |begins in an equal 
superposition |s) of all the vertices 

i 

\m) = I s ) = E 

Then it evolves by Schrodinger’s equation with Hamilto¬ 
nian 


H 0 = -jL - |w)H, 

where 7 is the probability amplitude per unit time of the 
particle jumping to an adjacent vertex (and is inversely 
proportional to mass), L is the graph Laplacian (an N- 
by-N matrix with 1 if two vertices are connected, 0 oth¬ 
erwise, and the negative of the degree of the vertex on 
the diagonal), and \w) is the “marked” vertex that we are 
searching for. The system evolves in a two-dimensional 
subspace spanned by {|to), |s)}. When 7 takes its critical 
value of 1/N, the energy eigenvectors are 

lfc ‘> = 7ivS <|s>±W) ' (1) 
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FIG. 1. (a) A strongly regular graph (Paley graph) with 

parameters (9,4,1,2). The “marked” vertex is colored red, 
vertices “one away” (adjacent) are colored blue, and vertices 
“two away” are colored white, (b) Success probability as a 
function of time for search on a strongly regular graph (Paley 
graph) with parameters (101,50,24,25). 


with corresponding energy eigenvalues 


E 0 ,i = —1 T 


1 

7N’ 


( 2 ) 


So the system evolves from |s) to |tu) in time tt\/N/2. 

In [4j, we showed that a randomly walking quantum 
particle can search for a marked vertex on a strongly reg¬ 
ular graph, an example of which is shown in Fig. [laj with 
the same asymptotic behavior as on the complete graph. 
An example of this is shown in Fig. [Tb] for N = 101, 
and it finds the marked vertex with probability near 1 
at time 7ry/lOT/2 « 15.786, as expected. So even though 
strongly regular graphs are not complete (and even lack 
global symmetry), they are “complete enough” for the 
search to primarily evolve in its two lowest energy eigen¬ 
states, which take the form of ([!]) with corresponding 
eigenvalues {2J, up to terms of order 1/y/N at the criti¬ 
cal 7. We call such graphs that evolve according to ([I]) 
and ([2]) with error terms that tend to zero for large N 
sufficiently complete. 

Strongly regular graphs are not the only sufficiently 
complete graphs—the hypercube is as well. The hy- 
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FIG. 2. (a) The 4-dimensional hypercube. The “marked” 
vertex is colored red, vertices “one away” (adjacent) are blue, 
vertices “two away” are yellow, vertices “three away” are 
green, and the vertex “four away” is white, (b) Success prob¬ 
ability as a function of time for search on the 10-dimensional 
hypercube, which has 2 10 = 1024 vertices. 


percube is even “less complete” than strongly regu¬ 
lar graphs—whereas search on strongly regular graphs 
evolve in a three-dimensional subspace, search on the n- 
dimensional hypercube evolves in a (n + 1 )-dimensional 
subspace, which grows with N = 2 n . An example of this 
in four dimensions is shown in Fig. |2a[ where vertices 
that evolve identically are the same color. Nonetheless, 
the hypercube is complete enough for search to behave 
like on the complete graph; search on it also primarily 
evolves in its two lowest eigenstates, which take the form 
of ([l]) up to terms of order 1/n at the critical 7 [3] [5j 6 ]. 
As an example, the evolution of the success probability 
for the 10-dimensional hypercube is shown in Fig. |2b| 
while N = 2 10 = 1024 is large enough for the runtime to 
be near 7 t\/1024/2 ss 50.265, it is not large enough for 
the second peak to be near 37 t\/ 1024/2 « 150.80, or for 
the success probability to be near 1 . 

As we detailed in 0 1 , we can speed up quantum 
search (ie., search on the complete graph) by evolving 
by the nonlinear Schrodinger equation 


But this improved runtime comes at the novel expense of 
increasing the time-measurement precision necessary to 
catch the sudden spike in success probability. Including 
time-measurement precision as a physical resource, it is 
better to choose g more modestly to opt for a less ag¬ 
gressive algorithm that runs in 0{N 1 ^ A ) time, which is 
a significant— but not unreasonable—improvement over 
the IV 1 / 2 scaling of Grover’s algorithm. The cubic-quintic 
and loglinear nonlinearities, on the other hand, achieve 
constant-runtime solutions without increasing the time- 
measurement precision, as shown in Figs. [3b| and|3c| 

This raises the question of whether computing with 
the nonlinear Schrodinger equation ([ 3 ]) speeds up search 
on sufficiently complete graphs in the same way as it did 
on the complete graph. In this paper, we show that it 
can, depending on the precise nonlinearity and graph. 
This indicates that our model of nonlinear continuous¬ 
time (quantum) computation can retain its power even 
when the underlying graph is not complete. 

II. LINEAR SEARCH 

We begin by introducing notation to describe linear 
search on sufficiently complete graphs. First we assume 
the system exactly evolves with eigenstates and eigenval¬ 
ues of the form 0 and 0, respectively, showing that 
it finds the marked vertex with probability 1 in time 
7 t\/N /2. Then we introduce error in the eigenstates and 
show how it propagates to the success probability and 
runtime. Since search on strongly regular graphs evolves 
in a three-dimensional subspace, and search on the hy¬ 
percube evolves in a (n + l)-dimensional subspace, search 
on a general sufficiently complete graph may evolve in an 
M-dimensional subspace. That is, the N vertices of the 
graph can be grouped together in M sets rrit of size \rrii\, 
where all vertices in a set evolve identically. Then the 
equal superpositions of identically evolving vertices 


= [ Hp- gf(\ip(r,t)\ 2 ) ]i/j{r,t), (3) 


H 


where / is a real-valued function. Many physically real¬ 
istic, effective nonlinearities take this form. For example, 
when f(p) = p, (|3j) is the Gross-Pitaevskii equation with 
a cubic nonlinearity that describes Bose-Einstein conden¬ 
sates under certain conditions; when f(p) = p — p 2 , the 
cubic-quintic nonlinearity describes nonlinear Kerr me¬ 
dia with defocusing corrections; and when f(p) = lnp, 
the loglinear nonlinearity describes Bose liquids mi- 
By choosing 7 at its time-varying critical value, the 
energy eigenvalues remain equal to ([l]), and H is equal 
to H 0 multiplied by a time-dependent term. This can be 
interpreted as a rescaling of time, so the nonlinear sys¬ 
tem evolves identically to the linear algorithm, but with 
a rescaled time. For example, with the cubic nonlinearity 
f(j>) = p , we can choose the nonlinearity coefficient g to 
achieve a constant-runtime solution, as shown in Fig. 3a 


| mi) 



E W 


form an orthonormal basis {|mo), |mi),..., |mM-i)} for 
the M-dimensional subspace. Without loss of generality, 
we pick the marked site to be |mo) = |u>), so |too| = 1. 
In the case of a strongly regular graph with parameters 
(N, k, A, p), the basis states of the subspace are 


|w 0 ) 

|mi) 

|m 2 ) 


\w) 

1 

Vk 


Y 1 *> 

(x,iv)££ 


1 

y/N-k-l 


Y i x >’ 

(x,w)g£ 


which correspond to the marked vertex, vertices adja¬ 
cent to the marked vertex, and vertices not adjacent to 
the marked vertex. For the n-dimensional hypercube, we 
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FIG. 3. Success probability as a function of time for search on the complete graph with N = 100 and N = 1000 vertices (black 
solid and red dashed curves, respectively) at the critical 7 with a: (a) cubic nonlinearity and g = N — 1, (b) cubic-quintic 
nonlinearity and g = N — 1, and (c) loglinear nonlinearity and g = y/N /logIV. 


first label each of the N = 2 n vertices with an n-bit string 
\z\...z n ). Without loss of generality, we choose the 
marked vertex to be the string of all zeros: | w) = |0 ... 0). 
Then the vertices “one away” are bit strings with a sin¬ 
gle one (ie., with Hamming weight 1), the vertices “two 
away” are bit strings with two ones (ie., with Hamming 
weight 2), and so forth. Taking the superposition of ver¬ 
tices equally far from the marked vertex, we get 

/ \— 1/2 

I TO fc)=(/J W...z n ). 

Then the set {| m/) : k = 0,1,..., n} is an orthonormal 
basis for the (n + l)-dimensional subspace. 

In this {|mo), |mi),..., |tom-i)} subspace, the state 
\ip(t)) °f the system can be written as a linear combina¬ 
tion of the basis states: 


M -1 

IV>(0} = C 00l TO i)' 

2=0 

We assume that the system evolves in its two lowest en¬ 
ergy eigenstates, having the form ([I]). Then the ampli¬ 
tudes are 


where = 1 when i = 0 and 0 otherwise is the Kro- 
necker delta. Then the amplitudes are 


Ci{t) = 


1 

Vn 


cos {vw) + isin {vw) 


y/N 


COS 


(wv) 


i = 0 
i / 0 


Squaring them, the probabilities are 

^ cos 2 (^)+ sin 2 (^v) 


m = 


N 

\™i\ 

N 


COS^ 


(vw) > 


i = 0 
i / 0 


This reveals that the success probability |c 0 | 2 reaches 1 
at time ny/N /2. 

Now, to introduce error, a sufficiently complete graph 
has eigenstates and eigenenergies of the form ([I]) and 0, 
respectively, both up to terms of some order e that tends 
to zero for large N: 


\M = -7= 


y/N 


y/2 V y/N± 1 


(|s) ± \w)) (1 + 0(e)), (4) 


E 0 ,i = -1 =F 


y/N 


(i + o(0)- 


( 5 ) 


d(t) = {mi\ip(t)) = {m/e lHt \s) 

= {mi\^o){^o\s)e~ lEot + (m i \ip 1 )(ip 1 \s)e~ lElt . 

Using 0. the inner products are 


(Oo.ilO 


y/2 y y/N ± 1 V y/N) 


1 I y/N VN± 1 

y/2 V y/N± 1 y/N 


{mi\4>o,i) 


J_ / Vn 

v^V Vn±i 


( a/KT 
v VN 



For strongly regular graphs, e = 1 /y/N, and for the n- 
dimensional hypercube, e = 1/n. Propagating these er¬ 
rors through the previous calculations, we get probabili¬ 
ties 


1 


N [y/N 
t 


t 


MOrH^cos 2 -/= (1 + 0(e)) 

(1 + 0 ( 0 ) 


- sm 


y/N 


and 


1 /.\i2 m 2 

]<Vo(OI = ^T C0S 


y/N 


(i + o(0) 


(i + o(0), 


(i + o(0). 


Thus the system evolves from |s) to \w) in time 7 ry/N /2 
for large N, as expected. 


and 
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III. NONLINEAR SEARCH 

For the nonlinear algorithm, we subtract from 
Hq an additional nonlinear “self-potential” V{t) = 
gf(\tf(r, f)| 2 ), where / is a real-valued function, so 
that the system evolves according to the nonlinear 
Schrodinger equation ^ . Then for positive g this speeds 
up the buildup of probability amplitude. In the compu¬ 
tational basis, the self-potential is 

N-l 

V{t) =g /(l<#}| 2 ) !*)(*!• (6) 

2=0 


phase), which is unobservable. We want to show that 
there exists a critical 7 that causes the nonlinear evolu¬ 
tion to follow the same path as the linear evolution for 
large N, but with rescaled time r. That is, we want to 
show that g(fi> 2 — / 1 ) can be dropped in comparison to 
1 + g(fo — fi) when we use the linear evolution, but with 
t —» r. Assume for the moment that we can do this. 
Then we have 

H = - 7 L - (1 + gf 0 - 5 / 1 ) |m 0 )<m 0 | 

for large N. Then the critical 7 for the nonlinear algo¬ 
rithm is 


Even with this nonlinearity, the system still 
evolves in the M-dimensional subspace spanned by 
{(too), |toi), ..., \rriM-i)}- In this M-dimensional sub¬ 
space, the self-potential ([b]) has off-diagonal elements 
equal to zero. Its diagonal terms are 


N -1 


(mi\V(t)\mi) = g^2 /(lO'IV’)l 2 ) {mi\j){j\mi) 

i=o 


jem; 


= 5 1 - 1 \mi\f 


= 9f Pr ■ 

\\m\J 

For ease of notation, let us define 

12 


I °i\ 
ImA 


fi = f Pi 

ImA 


Then in the M-dimensional subspace, the nonlinearity is 
(fo 0 ••• 0 \ 

0 h ••• i 


v(t) = g 


\o 


0 

0 fhi-ij 


In the M-dimensional subspace, the equation of motion 
is governed by 


M-l 

H = - 7 L - |m 0 )(m 0 | - g ^ fi\mi)(mi\ 

2=0 

M-l 

= -7 L ~ (1 +9fo) \m 0 ){mo\ - g ^ /i|mi)(m<| 

2=1 

= - 7 L - (1 + gf 0 - 5 / 1 ) |TO 0 )(m 0 | - gfil 

M-l 

- 9^2 (fi~ fi) I mi)(mi\. 

i—2 

The term proportional to the identity matrix can be 
dropped since it is a rescaling of energy (or an overall 


7c = 7 l (! + gfo - gfi ), 

where 7 l is the linear algorithm’s critical 7 . At 7 C , we 
have 


H = (I+ 5/0 - 5 / 1 ) (-7 lL - |m 0 )(m 0 |), 


which is the linear Hamiltonian at its critical 7 with a 
rescaled factor that depends on / 0 — f\. Thus the critical 
7 causes the nonlinear algorithm to evolve identically to 
the linear algorithm, but with rescaled time. For large 
N, the leading order behavior of fo — fi is 


fo - fl= cos 2 



+ sin 2 


/ 


iv cos 




which is the same for all sufficiently complete graphs that 
can be sped up by the nonlinearity, including complete 
graphs, and so we expect the nonlinearity to speed them 
up the same way. 

Now let us prove that such a critical 7 exists, namely 
that we can drop g(fi >2 — fi) compared to 1 + g(fo — / 1 ) 
with rescaled time t —> t in Eqs. Q and <©. beginning 
with the cubic nonlinearity f(p) = p. At later time, the 
sine term in fo clearly dominates both /)>2 and / 1 . To 
determine what happens at early time, we start with 


/i#o = jy cos 


Vn 


Cl + 0(c)) 


(1 + O(e)). 


Note that cos 2 (a(l + x)) ~ cos 2 (a) — 2acos(a) sin(a)a; for 
small x. Then we have 


~ N . 


- 2 


cos“ 


: COS 


Vn 


Vn ^ V/N 


sin 


Vn 


0(e) (1 + 0 (e)). 


For short time, cos (r/y/N) ~ 1 and sin (t/i/N) ~ 
t/\/N. Also absorbing —2 into the big-O, we get 


/m 0 




O(e) 


(1 + 0 (e)). 
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(a) (b) (c) 

FIG. 4. Success probability as a function of time for search on strongly regular graphs with parameters (IV, k, A, /r) = 
(509,254,126,127) and (4001,2000,999,1000) (black solid and red dashed curves, respectively) at the critical 7 with a: (a) cubic 
nonlinearity and g — N — 1, (b) cubic-quintic nonlinearity and g = N — 1, and (c) loglinear nonlinearity and g = VjV/ logIV. 



FIG. 5. Success probability as a function of time for search on the 10-dimensional and 11-dimensional hypercubes (black solid 
and red dashed curves, respectively) at the critical 7 with a: (a) cubic nonlinearity and g = N — 1 and (b) cubic-quintic 
nonlinearity and g = N — 1. Subfigure (c) is search on the 10-dimensional hypercube with g = 0 (linear), g = 0.01, g = 0.05, 
and g = 0.5 (black solid, red dashed, green dotted, and wildly varying blue solid curves, respectively). 


Multiplying out the terms, using cos (r/y/N) ss 1, and 
taking r = 0(1) we get, 

C0 ‘ 2 (tn) + 0{t) . ■ 

Then the terms we want to drop are, at short time, 



we have 


fi^O 



O(e) 


1 

N* 



+ 0(e) 


So the terms we want to drop are 


g(fi> 2 - A)« o (|). 

This can be dropped compared to 1 + g(fo — fi) when 
g = 0(jV) for both strongly regular graphs and the hy¬ 
percube, so the constant-runtime algorithm exists for 
both, as shown in Figs. |4a|and[5a| 

These arguments persist when f(p) = p q —the sine 
term in /o dominates fi^o at later time, and we have 
at early time 

cos2, (A) +0(e) .' 

Using this for the cubic-quintic nonlinearity f(p) = p—p 2 . 


g(fi>2 - h) = o (|), 

which can be dropped compared to 1 + g(fo — fi) when 
g = Q(N) for both strongly regular graphs and the hy¬ 
percube. So the constant-runtime algorithm exists for 
both, as shown in Figs. [4b] and Fig. |5b| Note that the 
peak in success probability is not as wide as for the com¬ 
plete graph in Fig.[3bJ this is expected because the width 
broadens as the success probability approaches 1. That 
is, as it approaches 1, / 0 goes to zero because of can¬ 
cellation between the cubic and quintic components, as 
does /1 because the failure probability goes to zero. Then 
the rescaling of time ceases, causing a broad peak. The 
error with which the success probability approaches 1 de¬ 
creases as N increases, so the peak is wider for larger N, 
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FIG. 6. Evolution of search with a loglinear nonlinearity and 
g = */N / log IV on a strongly regular graph with parame¬ 
ters (N, k, A, fj.) = (509,254,126,127). The black solid curve is 
|co(t)| 2 , the red dashed curve is |ci(f)| 2 , and the green dotted 
curve is |c 2 (t)| 2 ; they correspond to the marked vertex, ver¬ 
tices adjacent to the marked vertex, and vertices not adjacent 
to the marked vertex, respectively. 


as seen for strongly regular graphs in Fig. |4b| For the 
hypercube in Fig. |5b[ the success probability is too low, 
so we do not expect the peak to broaden at all. 

Now consider the loglinear nonlinearity f(p) = logp. 
Using the short time approximations from before and 
logo — log b = log (a/b), we have 


fi> 2 — fi 



Note that l/(a + x) ~ 1/a + (l/a 1 2 3 4 )a; for small x. Then 
we get 


fi> 2 ~ fi ~ log 



+ 0(e) 


x 




Multiplying this out and absorbing minus signs into the 
big-O’s, 


fi> 2 — /i 


log 




For short time, cos (r/y/N) ss 1, so 

fi >2 - fi ~ log (1 + O(e)). 

Since log(l + x) ~ x for small x, the terms we want to 
drop are 


g(fi >2 - fi) « O(ge). 


The constant-runtime algorithm requires g = 
&{VN/ log IV), so we can drop this compared to 
1 + g(fo — /i) for strongly regular graphs, as shown 
in Fig. [4c] The second “peak” is strange because of 
numerical error; the derivative of logs at x = 0 is 
nonzero, which makes the nonlinearity highly susceptible 
to noise, as shown in Fig. [6] where the evolution begins 
to vary wildly shortly after the first peak. For the 
hypercube, the constant-runtime algorithm does not 
exist. In fact, g must scale less than n = log IV, which is 
not significant enough for the loglinear nonlinearity to 
provide substantial speedup, as shown in Fig. [5c| 

So whether a sufficiently complete graph can be sped 
up by the nonlinear Schrodinger equation depends on the 
graph and the nonlinearity. Nonetheless, we have shown 
that even with a degree of noncompleteness, some nonlin¬ 
earities speed up search on certain sufficiently complete 
graphs in the same way as on the complete graph for 
large N. Thus completeness is not a requirement for fast 
nonlinear (quantum) search. 
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